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I. INTRODUCTION 



The experimental measurements and theoretical studies of the two body charmless 
hadronic B meson decays play an important role in the precision test of the standard 
mode (SM) and in searching for the new physics beyond the SM P|. For these charmless 
B meson decays, the dominant theoretical error comes from the large uncertainty in 
evaluating the hadronic matrix elements {MiM2\0i\B) where Mi and M2 are light final 
state mesons. The QCD factorization (QCDF) approach and the perturbative QCD 
(pQCD) factorization approach 0, 0, 13 are the popular methods being used to calculate 
the hadronic matrix element. Many charmless B meson decays, for example, have been 
calculated in the QCDF approach 12, iQl, ill, iSl] and/or in the pQCD approach [13, UlL 

SQIIl. 



In this paper, we would like to calculate the branching ratios and CP asymmetries for 
the four B ■nrj'^'^ decays by employing the low energy effective Hamiltonian and the 
pQCD approach. Besides the usual factorizable contributions, we here are able to evaluate 
the non-factorizable and the annihilation contributions to these decays. Theoretically, the 
four B wr^''> decays have been studied before in the naive or generalized factorization 
approach 0| as well as in the QCDF approach On the experimental side, the CP- 
averaged branching ratios and CP-violating asymmetries of S — > tt^t] and tt^t]' decays 
have been measured very recently 0, , and the world averages as given by the Heavy 
Flavor Averaging Group [l^ are 

Sr(S+ ^ vr+r;) = (4.3 ± 0.4) x 10"^ , (1) 

Br{B^ ^ vr+V) = (2.53l°|^) x 10"^ , (2) 

^cp(5^ ^ vr=^^) = -0.11 ±0.08 , (3) 

Acp{B^ ^ -n^r]') = 0.14 ±0.15. (4) 

For B — * 7r°?7, tt%' decays, only the experimental upper limits for the branching ratios are 
available now [19| 

Sr(S° n%) < 2.5 x 10"^ fir(5° ^ ttV) < 3.7 x 10^^ (5) 

For B nr]^''^ decays, the B meson is heavy, setting at rest and decaying into two light 
mesons (i.e. vr and r]^'^) with large momenta. Therefore the light final state mesons are 
moving very fast in the rest frame of B meson. In this case, the short distance hard process 
dominates the decay amplitude. We shall demonstrate that the soft final state interaction 
is not important for such decays, since there is not enough time for light mesons to 
exchange soft gluons. Therefore, it makes the pQCD reliable in calculating the B — > vrr/*^') 
decays. With the Sudakov resummation, we can include the leading double logarithms for 
all loop diagrams, in association with the soft contribution. Unlike the usual factorization 
approach, the hard part of the pQCD approach consists of six quarks rather than four. 
We thus call it six-quark operators or six-quark effective theory. Applying the six-quark 
effective theory to B meson decays, we need meson wave functions for the hadronization of 
quarks into mesons. All the collinear dynamics are included in the meson wave functions. 

This paper is organized as follows. In Sec. IHl we give a brief review for the PQCD fac- 
torization approach. In Sec. IIIIl we calculate analytically the related Feynman diagrams 
and present the various decay amplitudes for the studied decay modes. In Sec. I1V[ we 
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show the numerical results for the branching ratios and CP asymmetries oi B ^ ttt]^'^ de- 
cays and compare them with the measured values or the theoretical predictions in QCDF 
approach. The summary and some discussions are included in the final section. 



II. THEORETICAL FRAMEWORK 

The three scale PQCD factorization approach has been developed and applied in the 



lactorization approacn nas Dee: 
non-leptonic B meson decays B S S U) 1^ for some time. In this approach, 
the decay amplitude is separated into soft, hard, and harder dynamics characterized by 
different energy scales (t,mb, Mw)- It is conceptually written as the convolution, 

A{B ^ M1M2) ~ Jd%d%d% Tr [CmB{ki)^MAk2)^M,{h)H{kuk2,h,t)] , (6) 

where fcj's are momenta of light quarks included in each mesons, and Tr denotes the trace 
over Dirac and color indices. C (t) is the Wilson coefficient which results from the radiative 
corrections at short distance. In the above convolution, C{t) includes the harder dynamics 
at larger scale than Mb scale and describes the evolution of local 4-Fermi operators from 
mw (the W boson mass) down to t ~ 0{^/AMb) scale, where A = MB—rrib. The function 
H{ki, ^2, ^3,^) describes the four quark operator and the spectator quark connected by a 
hard gluon whose is in the order of AM^, and includes the 0{^/AMb) hard dynamics. 
Therefore, this hard part H can be perturbatively calculated. The function $m is the 
wave function which describes hadronization of the quark and anti-quark to the meson 
M. While the function H depends on the processes considered, the wave function is 
independent of the specific processes. Using the wave functions determined from other 
well measured processes, one can make quantitative predictions here. 

Since the b quark is rather heavy we consider the B meson at rest for simplicity. It is 
convenient to use light-cone coordinate {p'^,P~,Pt) to describe the meson's momenta, 

P^ = ^(P°±P'), and PT = {p\p^). (7) 

Using these coordinates the B meson and the two final state meson momenta can be 
written as 

Mr Mr Mr 

Pi = -|(1,1,0t), P2 = ^(1,0,Ot), P3 = ^(0,1,Ot), (8) 

respectively, here the light meson masses have been neglected. Putting the light (anti-) 
quark momenta in B, n and t] mesons as ki, k2, and k^, respectively, we can choose 

fci = (xiP+,0,kiT), ^2 = (x2P2+,0,k2T), k^ = (0, X^R^ Mt) ■ (9) 

Then, the integration over k^ , /c^, and k^ in eq.® will lead to 

A{B TCT]^'^) ~ Jdxidx2dxsbidbib2db2b3db3 

■Tr [CmBixi, 6i)$.(a;2, 62) V (^3, b^)H{x,, t)St{xi) e'^W] (10) 
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where h is the conjugate space coordinate of kiT, and t is the largest energy scale in func- 
tion H{xi,bi,t). The large logarithms {Inmw/t) coming from QCD radiative corrections 
to four quark operators are included in the Wilson coefficients C(t). The large double 
logarithms (In^ Xi) on the longitudinal direction are summed by the threshold resumma- 
tion and they lead to St{xi) which smears the end-point singularities on Xj. The last 
term, e"'^*^*-*, is the Sudakov form factor which suppresses the soft dynamics effectively 
[2ll |. Thus it makes the perturbative calculation of the hard part H applicable at inter- 
mediate scale, i.e., Mb scale. We will calculate analytically the function H{xi,bi,t) for 
B 7cr]^''> decays in the first order in expansion and give the convoluted amplitudes 
in next section. 



Wilson Coefficients 

For B TTT]^'^ decays, the related weak effective Hamiltonian H^ff can be written as 



n-/ - 
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i=3 



We specify below the operators in Tie// ioi h ^ d transition: 



0]' 


= daYLuf:) 




0^ 


= daY^'^a ■ 


U(3j^Lbf3 , 


03 


= d^YLba ■ 


Eg' Qp^i^Lq'^ . 


04 


= daYLbfs ■ 


Eg' Q'p^f^Lq'a , 
Eg' Q'pl^^RQa . 


05 


= d^YLba ■ 


Eg' Q'f3lt^R(l'f3 . 




= d^YLbp ■ 


07 


= ^daYLba 


■ Eg' eg'g^7/.^g^ , 

■ Eg' (^q'q.'piM'p > 


Og 


= \d^YLbp 


■ Eg' ^q'Q'pll^RQa 

■ Eg' (^q'Q.'plM'a 
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= ^daYLba 




= \d^YLbp 



(12) 



where a and (3 are the S'[/(3) color indices; L and R are the left- and right-handed 
projection operators with L = (1 — 75), R = (I + 75). The sum over q' runs over the 
quark fields that are active at the scale /i = O(mfe), i.e., {q'e{u,d, s,c,b}). The PQCD 
approach works well for the leading twist approximation and leading double logarithm 
summation. For the Wilson coefficients Cj(/i) {i = 1, . . . , 10), we will also use the leading 
order (LO) expressions, although the next-to-leading order (NLO) calculations already 
exist in the literature . This is the consistent way to cancel the explicit /i dependence 
in the theoretical formulae. 

For the renormalization group evolution of the Wilson coefficients from higher scale to 
lower scale, we use the formulae as given in Ref. directly. At the high myy scale, the 
leading order Wilson coefficients Ci{Mw) are simple and can be found easily in Ref.jl^. 
In PQCD approach, the scale t may be larger or smaller than the rrib scale. For the case 
of rrih < t < mwi we evaluate the Wilson coefficients at t scale using leading logarithm 
running equations, as given in Eq.(Cl) of Ref.[9]. In numerical calculations, we use 

Us = 47r/[/?i ln(t^/AQ^^ )] which is the leading order expression with AgJ,^ = 193MeV, 

derived from Ag^^ = 250MeV. Here /3i = (33 — 2n/)/12, with the appropriate number of 
active quarks Uf. Uf = 5 when scale t is larger than mi,. 
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At asiven energy scale t 
in Ref. are 



rrib = 4.8 GeV, the LO Wilson coefficients Ci{mb) as given 



-0.2703, C2 = 
0.0085, Cg = 
-0.0090, 



1.1188, C3 = 0.0126, C4 = -0.0270, 
-0.0326, C^ = 0.0011, Cg = 0.0004, 
= 0.0022. 



(13) 



If the scale t < mi,, then we evaluate the Wilson coefficients at t scale using the input of 
Eq. (jl3p and the formulae in Appendix D of Ref.j^ for the case of = 4. 



B. Wave Functions 



In the resummation procedures, the B meson is treated as a heavy - lig ht system. In 
general, the B meson light-cone matrix element can be decomposed as 22, 2^ 



(2vr)^ 

i 



iki -z 



{Q\hMdp{z)\B{pB)) 



V2K 



(ki) 



ifi — f - 



5b 



;ki) 



(14) 



13a 



where n = (1,0, Ot), and v = (0, 1, Ot) are the unit vectors pointing to the plus and 
minus directions, respectively. From the above equation, one can see that there are two 
Lorentz structures in the B meson distribution amplitudes. They obey to the following 
normalization conditions 



(2^ 



B 



d^ki 
(2^ 



^B 



(ki 



0. 



(15) 



In general, one should consider these two Lorentz structures in calculations of B meson 
decays. However, it can be argued that the contribution of (f)B is numerically small 
"2^ 2^, thus its contribution can be numerically neglected. In this approximation, we keep 



minimum number of input parameters for wave functions. Therefore, we only consider 
the contribution of Lorentz structure 



B 



V2K 



+"^B)750B(k 



1), 



(16) 



in our calculation. We use the same wave functions as in Refs.BEilillii. In the next 
section, we will see that the hard part is always independent of one of the kf and/or 
, if we make some approximations. The B meson wave function is then the function of 
variable k^ (or kf) and k^. 



(t)B{ki , k^ 



(17) 



The wave function for dd components in vr meson are given as 

1 



$^(P,a;,C) = 



^(x) + mo0^(x) + C"^o (# - ^ ■ '^)4>l{x)\ ■ 
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The wave function for dd components of r/'-'^ meson are given as 
1 



c 



Vdd 



Xx) + m^'^'^^Jx) + Cml^''{i)i - V ■ n)4>Ux) , (19) 



where P and x are the momentum and the momentum fraction of r^^j, respectively. We 
assumed here that the wave function of r]^^ is same as the vr wave function. The parameter 
C is either +1 or —1 depending on the assignment of the momentum fraction x. The ss 
component of the wave function can be similarly defined. 

The transverse momentum k-^ is usually conveniently converted to the h parameter by 
Fourier transformation. The initial conditions of leading twist (j)i{x), i = B,'n',ri,ri', are 
of non-perturbative origin, satisfying the normalization 



b = 0)dx = — -=fi 
2v6 



(20) 



with fi the meson decay constants. 



III. PERTURBATIVE CALCULATIONS 

In the previous section we have discussed the wave functions and Wilson coefficients of 
the amplitude in eq.®. In this section, we will calculate the hard part H{t). This part 
involves the four quark operators and the necessary hard gluon connecting the four quark 
operator and the spectator quark. We will show the whole amplitude for each diagram 
including wave functions. Similar to the B ^ ixp and B — > pri'^''> decays there 
are 8 type diagrams contributing to the B nrj^'^ decays, as illustrated in Figure 1. We 
first calculate the usual factorizable diagrams (a) and (b). Operators Oi, O2, O3, O4, O9, 
and Oio are {V — A){V — A) currents, the sum of their amplitudes is given as 

FeTT = SnCF'm'^ / dxidx^ / bidbib^db'i(f)B{xi,bi) 
Jo Jo 

X { [(1 + X3)0,(X3, 63) + (1 - 2s3)r,(0P(x3, 63) + 0t(x3, ^3))] 

■as{tl) he{xi, X3, bi, 63) exp[-Sab(tl)] 
+2r^(j)P{x3,b3)as(tl)he{x3,xi,b3,bi)exp[-Sabitl)]} , (21) 

where r^^ = 171^/7713', Cp = 4/3 is a color factor. The function he-, the scales t\ and the 
Sudakov factors Sab are displayed in Appendix|X| In the above equation, we do not include 
the Wilson coefficients of the corresponding operators, which are process dependent. They 
will be shown later in this section for different decay channels. 

The form factors of S to tt decay, F(fp'^(0), can thus be extracted from the expression 
in Eq. dH)), that is 

F^^iq' = 0) = Ff-(g2 = 0) = F,Jm|, (22) 

which is identical with that defined in Ref.0]. 

The operators O5, Og; O7, and Og have a structure of {V — A){y + A). In some decay 
channels, some of these operators contribute to the decay amplitude in a factorizable way. 
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FIG. 1: Diagrams contributing to the B vrr/^'^ decays (diagram (a) and (b) contribute to the 
B ^ TT form factor F^[^'^). 

Since only the axial- vector part of (V" + A) current contribute to the pseudo-sealer meson 
production, {7i\V - A\B){r]\V + A\0) = -{7i\V - A\B){r]\V - A\0), that is 

= -Fe. . (23) 

In some other cases, we need to do Fierz transformation for these operators to get right 
color structure for factorization to work. In this case, we get {S — P){S + P) operators 
from (y — A){V + A) ones. For these {S — P){S + P) operators. Fig. 1(a) and 1(b) give 

pi poo 

= — IGyrCirm^r^ / dxidx^ / 61^6163^63 0b(xi, 61) 
Jo Jo 

X { [0^(X3, 63) + r^((2 + X3)4)^{X3, 63) - X34)i{x3, 63)] 
■asitl)he{Xi, X3, 61, 63) exp[-S'afe(tg)] 

+ [xi(f)^{x3, 63) - 2{xi - l)r^(j)^{x3, 63)] 

■as{tl)he{x3,Xi,b3,bi)exp[-Sabitl)]} . (24) 
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For the non-factorizable diagrams 1(c) and 1(d), all three meson wave functions are 
involved. The integration of 63 can be performed using S function 5(63 — 62), leaving only 
integration of bi and 62- For the {V — A){V — A) operators, the result is 

Me-K = — ^ — nCpm^ / dxidx2dxz I 6id6i62C?&2 0s(a;i, &i)0r,(iC2, ^2) 
"J Jo Jo 

XX3 [2r^0^(a;3,62) - <^7r(.^3, ^2)] 

■as{tf)hf{xi, X2, X3, 61, 62) exp[-5'cd(t/)] ■ (25) 

For the {V — A){V + A) operators the formulae are different. Here we have two kinds 
of contributions from {V — A){V + A) operators. M^J and M^^ is for the {V — A){V + A) 
and {S — P){S + P) type operators respectively: 

Mf; = 0, (26) 
Mf,^ = -Me. . (27) 

The factorizable annihilation diagrams (g) and (h) involve only tt and r/*^'-* wave func- 
tions. There are also three kinds of decay amphtudes for these two diagrams. F^. is for 
{y — A){y — A) type operators, F^^ is for {V — A){y -\- A) type operators, while F^^ is 
for [S — P){S + P) type operators: 

p\ poo 

^ ^ ^^Cpm^ dX2dX3 / b2db2b3db3 {[X3(j)^{x3,b3)(j)r,{x2,b2) 

Jo Jo 

+2r^r.((x3 + 1)(/'^(X3, 63) + (xs - 1)0U^3, h))(t)^(x2, 62)] 
■ois{tl)ha{x2, X3, 62, bs) exp[-Sgh{tl)] 

- [X20.(X3,63)0^(X2,&2) 

+2r^r.(/)^(x3, 63)((a;2 + l)0j'(a;2, ^2) + (2^2 - 1)0* (2^2, ^2))] 
■a,{ti)haix3,X2,b3,b2)exp[-Sgh{tt)]} , (28) 



p1 poo 

F^J = 16TTCFm% / dx2dx3 / b2db2b3db3 
Jo Jo 

X { [X3r^{(f)^{x3, 63) - 0^(X3, 63))0^(X2, 62) + 2r^0.(x3, b3)(j)^{x2, 62)] 
XQ;s(t^)/la(x2, X3, 62, ^3) exp[-5'p^(t^)] 

+ [2r.0^(x3, 63)07,(2^2, 62) + X2rr,{(j)^{x2, 62) - 0ya:;2, &2))07r(a;3, ^3)] 
XQ;5(ie)^a(a:;3,2;2,&3,&2)exp[-,5gfe(i^)]} . (29) 

For the non-factorizable annihilation diagrams (e) and (f), again all three wave func- 
tions are involved. Here we have two kinds of contributions. M^., and describe 
the contributions from the {V - A){V - A), {V - A){V + A) and {S - P){S + P) type 
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operators, respectively, 

Ma-n = TrCpm'^ / dxidx2dx3 / 6ic/6i62'^&2 &l) 

3 Jo Jo 

X {- {X20vr(a;3, b2)(pr,{x2, 62) 

+r^r^ [{x2 + X3 + 2)0j(x2, 62) + (a;2 - a^s)^* (a;2, ^2)] 0^ (a^s, ^2) 
+ [(a;2 - X3)4>^{x3, 62) + (X2 + a;3 - 2)0* (x2, 62)] ^^(a^s, fe2)} 

■as{tj)hj{xi,X2, X3, bi, 62) exp[-S'e/(tj)] 
+ {a;30^(x3,62)0,y(a;2,&2) 

-r^r^ [0^(3:3, &2) [-(a;2 + X3)0^(x2, 62) + (a;2 - X3)4>l{x2, ^2)] 

-0t(a;3, ^2) - X2)4>^{X2, 62) + (a;2 + X3)0* (X2, 62)] ] } 

■asitj)hj{xi,X2,X3,bi,b2)exp[-Sef{tj)]} , (30) 

^OTT = — o — '^Cpm'^ / dxidx2dxs / 616/6162^^2 0_B(a;i, ^i) 

X { [(0:3 - 2)r^(f)^{x2, 62)(0f (x3, 62) + 0t(a;3, 62)) - (X2 - 2)r^0^(x3, 62) 
(0j'(x2, 62) + 0* (x2, 62))] ■ a,(tj)/ij(xi, X2, X3, 61, 62) exp[-S'e/(tJ)] 

- [a;3r^0^(x2, 62)(0^(X3, 62) + 0^(3^3, ^2)) 
-X2r^(f)^{x3, 62)(0^(X2, 62) + 0* (0:2, 62))] 

■as(t))/i)(a;i,X2,X3,6i,62) exp[-S'e/(t))]} , (31) 

_ _^y^^(j^^^ I (Ixidx2dx3 I 6i(i6i62(i62 0b(xi, 61) 
3 Jo Jo 

X { [a;307r(a;3, 62)0^(^2, 62) + r^r^(0^(a;3, 62)((2 + X2 + X3)0j'(x2, 62) 

-ix2 - X3)0i) + 0t(-(a;2 - a;3)0j + (-2 + X2 + ^3)0^))] 

■as{tj)^{xi, X2, X3, 61, 62) exp[-S'e/(tJ)] 

+ [-a;207r(a;3, 62)0^(^2, 62) - r^r^((x2 + ^3)0^(^2, 62) 

-(X2 - X3)0* (X2, 62)0^(0:3, 62) 

+ ((X3 - X2)0;^(X2, 62) - {X2 + X3)0^(X2, 62))0*^(X3, 62))] 

■a^(t))/i)(xi,X2,X3,6i,62)exp[-S'e/(t))]} . (32) 

In the above equations, we have assumed that xi « X2,X3. Since the light quark 
momentum fraction Xi in B meson is peaked at the small region, while quark momentum 
fraction X2 of r]^'^ is peaked around 0.5, this is not a bad approximation. The numerical 
results also show that this approximation makes very little difference in the final result. 
After using this approximation, all the diagrams are functions of /cf = xiuib/V^ of B 
meson only, independent of the variable of kf. Therefore the integration of eq. (|T7j) is 
performed safely. 

If we exchange the tt and r]^'^ in Fig. 1, the corresponding expressions of amplitudes for 
new diagrams will be similar with those as given in Eqs. ()21ll^ . since the vr and r]^'^ are all 



9 



pseudoscalar mesons and have the similar wave functions. The expressions of amphtudes 
for new diagrams can be obtained by the replacements 



For example, we find that 





0* - 








-F 


P^l - 














a-K ■ 



(33) 



where the form factors F^^{i), 



(34) 



and -F,^?,) describe the contributions induced by the 



iy - A)iy - A), iy - A)iy + A) and {S - P){S + P) operators, respectively. 

Before we put the things together to write down the decay amplitudes for the studied 
decay modes, we give a brief discussion about the rj—r]' mixing and the gluonic component 
of the r]' meson. 



The r] and rj' are neutral pseudoscalar (J^ = 0" 
mixtures of the SU{?>)f singlet rji and the octet rj^: 



mesons, and usually considered as 



cos 6p — sin Op 
sin On cos Or, 



(35) 



with 



^8 



—= [uu + dd — 2ss) 
V6 



[uu + dd + ss) 



(36) 



where Op is the mixing angle to be determined by various related experiments j28[. From 
previous studies, one obtains the mixing angle Op between —20° to —10°. One best fit 
result as given in Ref.|29] is -17° < 6'p < -10°. 

As shown in Eqs. fj35l36p . t] and t]' are generally considered as a linear combination of 



light quark pairs. But it should be noted that the r]' meson may has a gluonic component 
in order to interpret the anomalously large branching ratios of -B — > Kt]' and J/\l/ — 1]''^ 
"iil 3^ . In Refs.ji^, 3^, 3l|, the physical states i] and rj' were defined as 



\r] > 



\v' > 



Xr, 



uu + dd 



V2_ 
uu + dd 



V2 



+ Yff'lss > +Zr,i\gluonium >, 



(37) 



where X^(,) , and Z^f parameters describe the ratios of uu + dd, ss and gluonium 
{SU{?>)f singlet) componetnt of ri'^'\ respectively. In Ref. ^J, the author shows that the 
gluonic admixture in rj' can be as large as 26%, i.e. 



Zr,'/ (X,(,) + F^w + Zr,) < 0.26. 



(3^ 



10 



According to paper 0, a large SU(3) singlet contribution can help us to explain the 
large branching ratio for B — > Krj' decay, but also result in a large branching ratio for 
B K°ri decay, Br{B K^rj) ~ 7.0(13) x 10"^ for Op = -20°(-10°) as given in Table 
II of Ref. 3^ , which is clearly too large than currently available upper limits : 



Br{B K^rj) < 1.9 x 10~^ 



(39) 



Although a lot of studies have been done along this direction, but we currently still do 
not understand the anomalous gg — Vj coupling clearly, and do not know how to calculate 
reliably the contributions induced by the gluonic component of Vj meson. In this paper, 
we firstly assume that r]' does not have the gluonic component, and set the quark content 
of r] and r]' as described by Eqs. fj35l36|) . We will also discuss the effects of a non-zero 



gluonic admixture of r]' in next section. 

Combining the contributions from different diagrams, the total decay amplitude for 
B^ — i> 7r+?7 decay can be written as 



'7 5 2 1 1 1 

^'3 3 3 2 6 3 



— 6 ( C'a + 2^*4 — C's — gCe + -^Cj + gC's 



O 



Cs ) f:;F,{9,) 



-6 ( C, 



iuC2 — it 

1 



-c 



10 



-m2 it 



C3 

F2{e,) 
1 



2C4 — —Co 



-c 



10 



Fi[ 



2^6 + -Cs Fi(^p) + Ce - -Cs F^iO,) 



+ (M„^ + Me + Ma) [iuCl - it{C-i + Cg)] Fi(^^p) 

-(Mi;^+MfO 6(C^5 + Cr)-Fi(^^,) 
+ [F^^ + i^/'^ + i^r^) (-6) (Ce + Cs + ^C5 + \C,^ F, 



+Fe 



iuh;Ci + C2]-ith;C: 



Ca ~\~ —Co 

3 



C 



10 



Fi( 



f'F,{9,) 



-,^9 - -^Cio f^F2{6p) 



(40) 



where ^„ = V^^Kd, 6 = V^t^Vid, and 

Fiie^) 

F2{e,) 



— sin 6p + cos 

— sin Op — -\/2 cos ^. 



(41) 



are the mixing factors. The Wilson coefficients Cj should be_calculated at the appropriate 
scale t using equations as given in the Appendices of Ref. 
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Similarly, the decay amplitude for tt^t] can be written as 

1 



1 



it 



1 3 15 

-C3 — C4 — -C7 — -Cg + -Cg + Cio 



it ^3^3 + 2(^4 - 2C5 - -Ce - -C7 - -Cs + -Cg - gC'io^ ■ f^Fi{9p) 



'it ( C3 + 2^*4 — C5 ~ gC'e + + gC's ~ ~ 6^^° 



iuC2 — it[ + 2C4 — -Cg + 2^10 



-M 



-6 ( C4 - -Cio ) F^ie,) 



—it ( C'e + -C5 — -Cs — -C7 



1 3 

C3 + -Cg + 2^10 



+ {Ma^ + Ma + Me) iuC2 - it 



(42) 



The decay amplitudes for B 71^1]' and B tt^i]' can be obtained easily from 
Eqs. fl4()|l and ()42|1 by the following replacements 



rd rs ^ rd rs 

J Tj^ J Tj J r}' '> J rj' 

Flie,) - 

F2{e,) - 



f;(^,) = cos^, + ''''^^ 



F^iOp) = cos Op - V2 sin 9p. (43) 
Note that the possible gluonic component of t]' meson has been neglected here. 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



A. Input parameters and wave functions 

We use the following input parameters in the numerical calculations 



A. 



MS 



250MeV, 



m, 



l^^ = 1.4GeV, 



Mf 



5.2792GeV, M, 



w 



fn- 

fK 

80.41GeV. 



ISOMeV, 
= 2.4GeV, 



190MeV, 
= leOMeV, 



(44) 



For the CKM matrix elements, here we adopt the wolfenstein parametrization for the 
CKM matrix up to 0{X^), 



1 - 



CKM 



AX^{1 - p-ir]) 



A AX^{p - 17]) 
1 



(45) 



with the parameters A = 0.22, A = 0.853, p = 0.20 and r] = 0.33. 
For the B meson wave function, we adopt the model 



(l)B{x,b) = Nbx'^{1 — xYexp 



Ml 



(46) 



where iO\, is a free parameter and we take iOh = 0.4 ± 0.05 GeV in numerical calculations, 
and Nb = 91.745 is the normalization factor for uJt = 0.4. This is the same wave functions 
as m Refs.SIIilillii, which is a best fit for most of the measured hadronic B decays. 

For the light meson wave function, we neglect the b dependant part, which is not 
important in numerical analysis. We choose the wave function of vr meson (3^ : 



3 



(f)^{x) = -^f^x{l - X 



0.44^5/^ (2x 



0.25Cf^(2x 



2^6 



2V6 



1 + 0.43^2/^(20; - 1) + 0.09C]/^(2x - 1) 
1 - 2x) [1 + 0.55(10a;2 - lOx + 1)] . 



(47) 
(48) 
(49) 



The Gegenbauer polynomials are defined by 



d^^t) = i(3t2 - 1), C'j\t) = i(35t4 - 30t2 + 3), 

^(21t^-14t2 + l). 



he 



cl'\t) 



(50) 



For T] meson's wave function, 0^^^-, 0^^- and represent the axial vector, pseudoscalar 
and tensor components of the wave function respectively, for which we utilize the result 
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from the light-cone sum rule including twist-3 contribution: 

, . 3 



Vdd 



x] 



fxX{l 



+ar^4 [21(1 



x] 



1 

+8 



■fx 



1 



|l + a^| [5(1-2x^-1] 
14(l-2x)2 + l]|, 
[3(1 



2x) 



27 



11 



2a;)2 - 1] 
[35(1 -2x)^ 



30(1 - 2x)2 + 3] 



'ndd 



x] 



1 



2x) 



20 



:pl^_a1''^){Wx'^ - 10a; + 1) 



with 



,Vdd 



PVd 



0.44, 

mT,/m, 



.■ndd 



0.25, 



0.20, 



0.25, 



' 



773 = 0.015, CJ3 = -3.0. 



(51) 



(52) 



We assume that the wave function of uu is same as the wave function of dd. For the 
wave function of the ss components, we also use the same form as dd but with uiq^ and 
instead of m!^ and fx, respectively. For f^ and fy, we use the values as given in Ref. [3J| 
where isospin symmetry is assumed for f^ and S'f/(3) breaking effect is included for fy-. 



fx — f-K, fy — \ r^fic 



(53) 



These values are translated to the values in the two mixing angle method, which is 
often used in vacuum saturation approach as: 



/s = 169MeV, /i = 151MeV, 

08 = -25.9°(-18.9°), 6i = -7.r(-0.1°) 

where the pseudoscalar mixing angle 9p is taken as —17° (—10°) |2 
= Vdd{uu),Vss) are defined as: 

i2 



^dd{uu) 





ml 



rrin 



2Ml 



m: 



(54) 

The parameters itIq 
(55) 



(m„ + mrf)' """^ (2m J 

We include full expression of twist— 3 wave functions for light mesons. The twist— 3 
wave functions are also adopted from QCD sum rule calculations ji^. We will see later 
that this set of parameters will give good results for B — ^ tit]^'^ decays. Using the above 
chosen wave functions and the central values of relevant input parameters, we find the 
numerical values of the corresponding form factors at zero momentum transfer: 



(56) 



pB- 




= 0) 


= 0.30, 


pB- 




= 0) 


= 0.15, 


pB- 




= 0) 


= 0.14. 
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These values agree well with those as given in Refs. 



B. Branching ratios 



For B — >• TTT]^'^ decays, the decay amplitudes in Eqs. (jiUI) and can be rewritten as 



where 







p 


V*hVud 




T 



(57) 



(5^ 



is the ratio of penguin to tree contributions, a = arg — t/''v* is the weak phase (one of 

the three CKM angles), and S is the relative strong phase between tree (T) and penguin (P) 
diagrams. The ratio z and the strong phase 6 can be calculated in the pQCD approach. 
One can leave the CKM angle a as a free parameter and explore the CP asymmetry 
parameter dependence on it. 

For B — s> TT^f] decay, for example, one can find "T" and "P" terms by comparing the 
decay amplitude as defined in Eq. (pOj) with that in Eq. (jSTj), 



FA 



Ci + -C2 ) + M,^C2 



+Fe -Ci + C2 + (M, + Me + Ma.) C, 



(59) 



+ ( C3 + - C's - + -^C^ + -Cg - -Cg - -C'lo ) f^F2{6p) 



P2 



— -\- Ca — —C' 

3 6 



+m5 



~ ( C*3 + 2C4 — -Cg + 2^10 ) Fl{^p/ 



10 J^2K^p) 



+ {Ma. + M, + MJ(C3 + Cg)Fi(^p) + {m2 + Mf^) {C, + Cr) ■ F,( 
HF^^ + F^- + Fr^)Fi(^,) (^Ce + C, + \c, + ^C^) 



+Fe ( ^Ca + C4 + ^Cg + Clo ) 



(60) 



Similarly, one can obtain the expressions of the corresponding tree and penguin terms for 
the remaining three decays. 
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Using the "T" and "P" terms as given in Eqs. and (jfiOj) . it is easy to calculate the 
ratio z and the strong phase b for the decay in study. For p^r] and p^rj' decays, 

we find numerically that 

z{7r+ri) = 0.33, ^(vr+r]) = -136°, (61) 
^(vr+r/') = 0.25, (5(7r+V) = -130°. (62) 

The main errors of the ratio z and the strong phase 6 are induced by the uncertainty 
of ujf, = 0.4 ± 0.05 GeV and ttiq = 1.4 ± 0.1 GeV and small in magnitude. The reason 
is that the errors induced by the uncertainties of input parameters are largely cancelled 
in the ratio. We therefore use the central values of z and 6 in the following numerical 
calculations, unless explicitly stated otherwise. 

From Eq. ()57|) . it is easy to write the decay amplitude for the corresponding charge 
conjugated decay mode 

M = V^,V:,T - Vt,V;,P = VMT [I + ze^(-"+')] . (63) 

Therefore the CP-averaged branching ratio for nr]^'^ is 

Br = {\M\^ + \M\^)/2= KbV*aT\^ [l + 2z cos a cos 6 + z^] , (64) 

where the ratio z and the strong phase 6 have been defined in Eqs. ()57p and (j^Hj) . It is easy 
to see that the CP-averaged branching ratio is a function of cos a. This gives a potential 
method to determine the CKM angle a by measuring only the CP-averaged branching 
ratios with PQCD calculations. 

Using the wave functions and the input parameters as specified in previous sections, it 
is straightforward to calculate the branching ratios for the four considered decays. The 
theoretical predictions in the PQCD approach for the branching ratios of the decays under 



consideration are the following 

Br{ B^ TT+r^) = [4.1+J|(c^,)l°:^(m^)^°i(a)] x 10-^ (65) 

Br{ B+ n+ri') = [ZAtoHcOb) ± 0.2(m^) ± 0.3(a)] x 10-^ (66) 

Br{ B' A) = [0.23l°:°^(o.,)+°:°^(m^) ± 0.05(a)] x 10-^ (67) 

Br{ 5° ^ ttV) = [0.19 ± 0.02(cjfe) ± 0.03{m^)tom{a)] x 10~\ (68) 

for 9p = -17°, and 

Br{ B+ n+ri) = [S.StHitOb) ± 0.3(mJ) ± 0.4(a)] x 10"^ (69) 

Br{ 5+ ^ TT+r/') = [3.2tl-^j{ub) ± 0.3(mJ) ± 0.4(a)] x 10'^, (70) 

Br{ 5° ^ A) = [0.17 ± 0.02{ujb) ± 0.02{m^)toM{a)] ) x 10-^ (71) 

Br{ 5° ^ ttV) = [0.28tom{^b) ± 0.04(mJ) ± 0.05(a)] x 10~^ (72) 



for 6p = —10°. The main errors are induced by the uncertainties of ujb = 0.4 ± 0.04 GeV, 
= 1.4 ± 0.1 GeV and a = 100° ± 20°, respectively. 

The PQCD predictions for the branching ratios of considered decays agree well with 
the measured values or the upper limits as shown in Eqs.(IT]|HI). Furthermore, the pQCD 
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FIG. 2: The a dependence of the branching ratios (in unit of 10 ^) of vr+r/^'^ decays 

for ml = 1.4 GeV, Op = -17°, = 0.36 GeV (dotted curve), 0.40 GeV (sohd curve) and 0.44 
GeV(dashed curve). The gray band shows the experimental data. 



predictions also agree well_with the theoretical predictions in the QCDF approach, for 
example, as given in Ref. 



Br{ 5+ - 
Br{ fi+ - 
Br{ 5° 
Br{ B^ - 



+ i\ 
^ V) 

.A) 



1-6 



1-6 



(3.lli;?) X 10- 
(0.28l°:i) X 10-^ 

(0.17^°:??) X 10"^ 



(73) 



where the individual errors as given in Ref. p] have been added in quadrature. 

It is worth stressing that the theoretical predictions in the PQCD approach have rel- 
atively large theoretical errors induced by the still large uncertainties of many input 
parameters, such as Ub, mg and 9p. In our analysis, we consider the constraints on these 
parameters from analysis of other well measured decay channels. For example, the con- 
straint 1.1 GeV < ttIq < 1.9GeV was obtained from the phenomenological studies for 
B ^ nn decays while the constraint of a ~ 100° ± 20° were obtained by direct mea- 
surements or from the global fit 0, |3^. From numerical calculations, we get to know 
that the main errors come from the uncertainty of Uh, mg, a and 9p. 

In Figs. 121 and El we present, respectively, the pQCD predictions of the branching ratios 



of 5 



TC 



'rj and tt^^t]' decays for 6p = 10° 



0.4 ±0.04 GeV, 



rrin 



1.4 ±0.1 GeV and 



a = [0°, 180°]. Fig. m shows the a-dependence of the pQCD predictions of the branching 
ratios of B ^ vt^t^W decays for Op = 10°, Ub = 0.4 GeV, = 1.4 GeV and a = [0°, 180°]. 

From the numerical results and the figures we observe that the pQCD predictions 
are sensitive to the variations of Ub and mg. The parameter ttIq originates from the 
chiral perturbation theory and have a value near 1 GeV. The ttiq parameter characterizes 
the relative size of twist 3 contribution to twist 2 contribution. Because of the chiral 
enhancement of ttiq, the twist 3 contribution become comparable in size with the twist 2 
contribution. The branching ratios of Br(B nr]^'^) are also sensitive to the parameter 
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FIG. 3: The a dependence of the branching ratios (in unit of 10 ^) of — > ■n^rj-''^ decays for 
= 0.4 GeV, 6'p = -17°, = 1.3 GeV (dotted curve), 1.4 GeV (soHd curve) and 1.5 GeV 
(dashed curve). The gray band shows the data. 
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FIG. 4: The a dependence of the branching ratios (in unit of 10 ®) of vr'^r; (sohd curve) and 
7r°?7' (dotted curve) decays for = 1.4 GeV, = -17°, = 0.40 GeV . 

rriQ, but not as strong as the dependence. 



C. CP- violating asymmetries 

Now we turn to the evaluations of the CP-violating asymmetries of i? ^ Tir]^''^ decays 
in pQCD approach. For tt^t] and ti^t]' decays, the direct CP-violating 

asymmetries Acp can be defined as: 



2z sin a sin (5 



|y\/l|2 + |_/\^|2 l + 22cosacos5 + z2' 



(74) 
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where the ratio z and the strong phase 8 have been defined in previous subsection and 
are calculable in pQCD approach. 

It is easy to calculate the CP-violating asymmetries with z and b. In Fig. we show 
the a— dependence of the direct CP-violating asymmetries .4^*p for n^ri (the solid 

curve) and — > vr^r]' (the dotted curve) decay, respectively. From Fig. |31 one can see 
that the CP- violating asymmetries J^Qp{B^ n^i]^'^) are large in magnitude, about 35% 
for a ~ 100°. So large CP- violating asymmetry plus large (~ 10~^) branching ratios are 
measurable in current B factory experiments. 

The pQCD predictions for A'^p and the major theoretical errors for B"^ -k^t]^'^ 
decays are 

A'^^p{B^ TX^V) = (-37l^(a) ± 4(u;,)i;(m-)) x 10-^ (75) 
A'SUB^ - nS) = (-33^^(a)l^(a;,)l°(m^)) x 10-^ (76) 

where the dominant errors come from the variations of Ub = 0.4±0.04 GeV, m'^ = 1.4±0.1 
GeV and a = 100° ± 20°. 

By comparing the above numerical results with those measured values as given in 
Eqs.dSD and \ we find that 

1. The pQCD predictions for the direct CP- violating asymmetry for both B ^ nr] and 
TTT]' decays are large in magnitude and have a moderate theoretical error because of 
the cancelation in the ratios. 

2. For B^ TT^i] decay, the pQCD prediction for A^p has the same sign with the 
measured value and also consistent with it within 2a errors. For tt^i]' decay, 
however, the pQCD prediction for A^p has the opposite sign with the measured 
value. 

3. For B^ n^r]' decay, although there exist a clear difference between the pQCD 
prediction of A^p and the data, but it is too early to draw any reliable information 
from such difference because of the still large theoretical and experimental errors. 
More theoretical studies (for example, calculation of next-to-leading order contribu- 
tions [3^) and more accurate measurements are needed to clarify this discrepancy. 

In Ref. 0, by using the "default values" of input parameters, the authors presented 
their predictions for Aqp{B -K^r]^'^) in the QCDF approach 

A^^p{B^ - ^^V) = {-l^.9tti tli X 10-, (77) 
A'SUB^ n^v') = {-S.Qtli tZ' trrtlVd x 10^, (78) 

where the first error comes from the variation of the CKM parameters, the second error 
refers to the variation of /i ~ m;,, quark masses, decay constants, form factors, and the 
mixing angle 6p. The third error corresponds to the uncertainty due to the Gegenbauer 
moments in the expansion of the LCDAs. The last but largest error is induced by the 



^ We used the same convention as the BaBar and BeUe Cohaboration [13, ll^ to define the CP-violating 
asymmetry for B ttt]^'^ decays. 
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FIG. 5: The direct CP asymmetries (in percentage) of 
(dotted curve) as a function of CKM angle a. 



Tr~^ rj (solid curve) and B 



uncertainty of the unknown annihilation contributions. In fact, these numbers are similar 
to those generalized factorization approach since the mechanism of strong phase is 
the same for these two approaches. 

By using the "set 5*4" input parameters jsl, the central values of the QCDF prediction 
can become positive simultaneously 



5.6 X 10"^ 
11.1 X 10" 



(79) 
(80) 



From the numerical results as given in Eqs. 1)7711^0)) . one can see that the dominant 
theoretical error from the annihilation contributions in QCD factorization approach are 
too large to make any meaningful comparisons between the theoretical predictions and 
the data for ^^^^,(5 ^ ttt]^'^). The reason is that the annihilation contributions play a 
key role for producing the strong phase of the two-body charmless B meson decays, but 
unfortunately they are incalculable in QCD factorization approach. 

We now study the CP-violating asymmetries for 
decay modes, the effects of B^ — 5° mixing should be considered. For 5° meson decays 
we know that AF/Am^ <^ 1 and AF/F <^ 1. The CP- violating asymmetry of B^{B^) — 
7^0^(0 (^gcay is time dependent and can be defined as 



7j-o^(') decays. For these neutral 



Acp = 



F SO(At)^/cP -r(SO(At)^/< 



CP) 



F (^i?o(At) ^ fcp) + F (i?o(At) ^ fcp) 
= A'^jp cos(AmAt) + A"^^ sin(AmAt), 

where Am is the mass difference between the two B^ mass eigenstates. At 



tcp " t 



is the time difference between the tagged 5° {B ) and the accompanying B (5°) with 
opposite b flavor decaying to the final CP-eigenstate fcp at the time tcp- The direct and 



tag 
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mixing induced CP-violating asymmetries and A^p can be written as 



Adir _ l^Cpf - 1 .mix _ '^Im{Xcp) , . 

•^'"p ~ 1 I I \ 12' CP — 1 I I \ F' y^-^i 
1 + \acp\ l + \Xcp\ 



where the CP- violating parameter Xcp is 

Vt,V;,{7rOr^(')\H,ff\BO) 1 + ze^(^+") " ^'^''^ 

Here the ratio z and the strong phase 6 have been defined previously. In PQCD approach, 
since both z and 6 are calculable, it is easy to find the numerical values of A'^p and A^p 
for the considered decay processes. 

In Figs, inland [71 we show the a— dependence of the direct and the mixing-induced CP- 
violating asymmetry and A^p for n^r] (solid curve) and ti^t]' (dotted 
curve) decays, respectively. 

The pQCD predictions for the direct and mixing induced CP- violating asymmetries of 
B^ Ti'^r]^'^ decays are 

A'SUB' - vr%) = {-A2tlM&uj,)tl{m-,)) x 10-^ (84) 

A'SUB' - vr^O = {-36tl\a)tl{uj,)tl{m-,)) x 10-^ (85) 

^^i?(i?° ^ 7r%) = {67tlia)tliuj,)tlirn-,)) x 10-^ (86) 

(S° ^ vr^O = {Q7tl{a)tt{u^,)tl{m-,)) x 10-^ (87) 

where the dominant errors come from the variations of uJh = 0.4±0.05 GeV, itlq = 1.4±0.1 
GeV and a = 100° ± 20°. 

As a comparison, we present the QCDF predictions for Aqp{B^ n^r]^'^) directly 
quoted from Ref. ||8|] 

A'SUB' - ir%) = {-17.9tli tlf/h't!^l.) X 10-, (88) 
A'SUB' - A') = {-I9.2tli til m nii) X 10-^ (89) 



where the "default values" of the input parameters have been used |8|, and the error 
sources are the same as those for the numerical results in Eqs. (j77p and (j78p . Currently, 
no relevant experimental measurements for the CP-violating asymmetries of i? tc'^t]^'^ 
decays are available. For the direct CP-violating asymmetries of 5" n'^r]^'^ decays, 
the theoretical predictions in pQCD and QCDF approach have the same sign, but the 
theoretical errors are clearly too large to make a meaningful comparison. One has to wait 
for the improvements in both the experimental measurements and the calculation of high 
order contributions. 

If we integrate the time variable t, we will get the total CP asymmetry for B'^ ti'^t]^'^ 
decays, 

AcP = Y^.A'cP + YT^2^CP^ (90) 

where x = Am/T = 0.771 for the B^ — B^ mixing fis']. In FiglHl we show the a-dependence 
of the total CP asymmetry Acp for B^ — > ti^t] (solid curve) and 5° — > tt^t]' (solid curve) 
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FIG. 6: The direct CP asymmetry A^'p (in percentage) of vr^r/ (solid curve) and ir^r]' 
(dotted curve) as a function of CKM angle a. 
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FIG. 7: The mixing induced CP asymmetry A^'p (in percentage) of B^ 
and B^ ir^rj' (dotted curve) as a function of CKM angle a . 



TT^T] (solid curve) 



decay, respectively. Numerically, we found 



Atot 



.cp{B' - A) = (+3.7t?i(«)l|^(^,)l?:^(m^)) X 10-^ 
Afp{B'^.%') = (+7.7i^:?(a)l°i(c.,)l^:?(m^)) X 10^1 



(91) 
(92) 



D. Effects of possible gluonic component of r]' 

Up to now, we have not considered the possible contributions to the branching ratios 
and CP-violating asymmetries of i? — *• ttt]' decays induced by the possible gluonic com- 
ponent of Tj' j^]. When Z^' 7^ 0, a decay amplitude M.' will be produced by the 
gluonic component of ry'. Such decay amplitude may construct or destruct with the ones 
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FIG. 8: The total CP asymmetry (in percentage) of — > -K^rj (solid curve) and tt^t]' 
(dotted curve) as a function of CKM angle a . 

from the qq {q = u, d, s) components of rj', the branching ratios of the decays in question 
may be increased or decreased accordingly. 

Unfortunately, we currently do not know how to calculate this kind of contributions 
reliably. But we can treat it as an theoretical uncertainty. For \M' /M{qq)\ ~ 0.1 — 0.2, 
for example, the resulted uncertainty for the branching ratios as given in Eqs. ()(i(i|l and 
(|68p will be around twenty to thirty percent. 

From Eq. ()2I66|) . one can see that the theoretical result of Br{B^ n^r]') in the 
PQCD approach agree very well with the measured values within one standard deviation. 
Furthermore, the pQCD predictions for the branching ratios of B ^ pr]^'^ decays also show 
very good agreement with the data . We therefore believe that the gluonic admixture 
of 1]' should be small, and most possibly not as important as expected before. 

As for the CP-violating asymmetries of -B ^ ttt]' decays, the possible contributions of 
the gluonic components of the t]' meson are largely cancelled in the ratio. 

V. SUMMARY 

In this paper, we calculate the branching ratios and CP-violating asymmetries of B^ — > 
TT^r], B^ TT^T]', B^ —>■ TT^T], and B^ —>■ -n^Vj decays in the PQCD factorization approach. 

Besides the usual factorizable diagrams, the non-factorizable and annihilation diagrams 
are also calculated analytically. Although the non-factorizable and annihilation contribu- 
tions are sub-leading for the branching ratios of the considered decays, but they are not 
negligible. Furthermore these diagrams provide the necessary strong phase required by a 
non-zero CP-violating asymmetry for the considered decays. 

From our calculations and phenomenological analysis, we found the following results: 

• From analytical calculations, the form factors for B ^ rj, B ^ t]' and B ^ n tran- 
sitions can be extracted. The PQCD results for these form factors are F^{^'^{0) = 

0.30, Fo^r''(0) = 0.15 and F^^^''' (0) = 0.14. 
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For the CP-averaged branching ratios of the four considered decay modes, the the- 
oretical predictions in PQCD approach are 

Br{B+ TT+r/) = (4.1+^;^) x 10-^ 

Br{B+ ^ A^'^) = (2.4l[J:^) X 10-^ 

Br{B°^n%) = (0.23 ± 0.08) x 10"^ 

Br{B+ ^ Ti^ri^'^) = (0.19 ± 0.05) x 10-^ (93) 



where the various errors as specified in Eqs. ()65ll68p have been added in quadrature. 
Although the theoretical uncertainties are still large (can reach 40%), the leading 
PQCD predictions agree very well with the measured values or currently available 
experimental upper limits, and are also well consistent with the results obtained by 
employing the QCD factorization approach. 

For the CP-violating asymmetries, the theoretical predictions in PQCD approach 
are 

A'SUB^ - ^^V) = (-371?) X 10-^ 
A'SUB^ - rr^v') = (-33^^) x 10-^ 
A'SUB' - A) = (-421}°) X 10-^ 
/l™"(5° ^ A) = (+67l?i) X 10"^ 
A'SUB' - A') = (-361}^) X 10-^ 

A^i? (5° ^ A') = (+67l?i) X 10-2, (94) 

where the various errors as specified in Eqs. ()75I76|) and ()84ll87p have been added 
in quadrature. The uncertainties are around twenty to forty percent. And finally, 
the time integrated CP asymmetry for the neutral decay modes are 

A'SUB' - 7r%) = {+3.7lli) X 10-^ 

A'Sp{B' ^ vr%') = (+7.712-7) X 10-2. (95) 

For B^ TT^i] decay, the pQCD prediction for A^p has the same sign with the 
measured value and consistent with the data within two standard deviations. For 
5^ — > n'^rj' decay, however, the pQCD prediction for A'^^p has an opposite sign 
with the measured value. Great improvements in both the theoretical calculations 
and experimental measurements are needed to clarify this discrepancy. 

From the very good agreement of the pQCD predictions for the CP-averaged branch- 
ing ratios Br{B^ iT~^ri^'^) and Br{B^ p'^rj^'^) 14] with the measured values, 
we believe that the gluonic admixture of 7]' should be small, and most possibly not 
as important as expected before. 
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APPENDIX A: RELATED FUNCTIONS 

We show here the function /ij's, coming from the Fourier transformations of H^'^\ 

heixi,X3, bi, 63) = Kq {y/x^niBbi) [9{bi - b3)Ko {y/x'^mBbi) Iq {y/xiuiBb^) 

+9{bs - bi)Ko (y/x^niBbs) Iq {^/x^mBbl)] St{x3), (Al) 

haix2, X3, 62, 63) = Ko {iy/x^niBb^) [d{b3 - b2)Ko {iy/x^rriBbs) Iq («v^"^i?^2) 

+e{b2 - 63)^0 iiVx^mBb2) lo (^yis^^s^s)] Stixs), (A2) 



/l/(Xi,X2,X3,6i,62) = \9{b2 - bl)lo{MB^/x^ibl)Ko{MB^/x^b2) 

Ko(MbF(i)6i), for F^) > 0^ 
+ {bi - 62) M ^H«(M^^7f^ b,), for Ff,) < 1 ' ^^^^ 



h^f{xi, X2, X3, bi, 62) = \ 0{bi - b2)KQ{i^X2X3biMB)h{i^Jx2X2,b2MB) + (61 ^ 62' 



2 '^o\^xi + X2 + X3 - X1X3 - X2a;3 61 M^), (A4) 



hj{xi,X2,X3,bi,b2) = ^(^1 ~ b2)Ko{i^/x^blMB)lo{i^/x^b2MB) 

>| / Ko(MbF(2)6i), for Ff2)>0\ 
+ ib. - 62)) • ^^H«(M,^ 6,), for Fg^ < j ' ^^5) 

where Jq is the Bessel function and Kq, Iq are modified Bessel functions KQ{—ix) = 
— {'7i/2)Yq{x) + i{TT/2)Jo{x), and F(j)'s are defined by 

^(1) = i^l - ^2)X3 , (A6) 

-^(2) = i^i - ^2)^3 ■ (AT) 
The threshold resummation form factor St{xi) is adopted from Ref . |2f 

^y7^l [l + C) 

where the parameter c = 0.3. This function is normalized to unity. 
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The Sudakov factors used in the text are defined as 

Sab{t) = s [xirriB/ V2, 61 j + s [x^niB/ v^, 63) + ^ (^(1 - x^^itib/ v^, 63) 
1 



-ln(6iA) -ln(63A) 
s {xiItlbI V2, 61 j + s {x2mB/ V2, 62 j + s ^(1 

1 



(A9) 



X2 yrriB 



/V2,h 



Pi 



21nim + ln.l-(^/^) 



ln(6iA) ' '"-ln(62A)_ 
s (^xitjib/ V2, 61 j + s {x2mB/ V^, 62 j + S ^(1 - X2)mB/ V^, 62 j 

+s (^X3mB/ v^, 62) + s (^(1 - x^)mB/ V2, ^2) 

A ['''-ln(6iA)^^'''-ln(62A)_ ' 

s {x2mB/ 61 j + s (^XsTJIb/ V2, 62) + s ^(1 - X2)mB/ V2, h 



(AlO) 



+ s 



X3)mB/V2,b2) - 



-ln(6iA) -ln(62A) 



where the function b) are defined in the Appendix A of Ref.j9|. The scale t, 
above equations are chosen as 



(All) 

(A12) 
's in the 



tl 


= max{y/x^mB, l/bi, I/63) , 


(A13) 


tl 


= max{y/x^mB, l/bi, I/63) , 


(AM) 




= max{y/x^mB, I/62, 1/&3) , 


(A15) 




= max(v/a^mB, I/62, 1/^3) , 


(A16) 


tf 


= max(v/xiX3mB, a/ (xi - X2)x^mB, I/&1, 1/&2) , 


(A17) 




= max(Vxi + X2 + X3 - X1X3 - X2X3mB, y/X2X3mB, I/&1, 1/&2) , 


(A18) 




= max(i/a;2a;3mB, I/61, 1/62) • 


(A19) 
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